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We present a computational study of individual and pairs of substitutional Mn impurities on the 
(110) surface of GaAs samples based on density functional theory. We focus on the anisotropy prop- 
erties of these magnetic centers and their dependence on on-site correlations, spin-orbit interaction 
and surface-induced symmetry-breaking effects. For a Mn impurity on the surface, the associated 
acceptor-hole wavefunction tends to be more localized around the Mn than for an impurity in bulk 
GaAs. The magnetic anisotropy energy for isolated Mn impurities is of the order of 1 meV, and 
can be related to the anisotropy of the orbital magnetic moment of the Mn acceptor hole. Typically 
Mn pairs have their spin magnetic moments parallel aligned, with an exchange energy that strongly 
depends on the pair orientation on the surface. The spin magnetic moment and exchange energies 
for these magnetic entities are not significantly modified by the spin-orbit interaction, but are more 
sensitive to on-site correlations. Correlations in general reduce the magnetic anisotropy for most of 
the ferromagnetic Mn pairs. 



I. INTRODUCTION 

Since the discovery of ferromagnetic order in Mn- 
doped GaAs with a Curie temperature above of 
100 Kfi research in dilute magnetic semiconductors 
(DMS) has developed into an important branch of 
material science. Early work, besides being aimed 
at understanding the physics of (Ga, Mn)As was 
strongly focused on the goal of achieving room- 
temperature ferromagnetism in this prototype DMS. 
Although this goal seems now of difficult realization, 
(Ga, Mn)As and DMSs in general still attract a lot of 
attention both for fundamental science and applica- 
tions (e.g. in spintronics)j2^ From the point of view 
of theory, intense effort based on both phenomeno- 
logical models^— and ab-initio calculations 1 ^ - — has 
lead to undoubted progress in understanding the ori- 
gin and the properties of ferromagnetism in (Ga, 
Mn)As. Yet, some fundamental aspects related to 
the microscopic mechanism remain controversial and 
are still strongly debated*^ - — 

On the experimental front, recent scanning tun- 
neling microscopy (STM) studies in (Ga, Mn) As 
have been able to visualize the spatial structure 
of the acceptor wavefunctions at doping concentra- 
tions near the metal-insulator transition and shown 
that it has a multifractal character, possibly indicat- 
ing carrier correlations^ Similar STM experiments 
have been carried in the last few years at very low 
doping concentrations^ - — probing specifically the 
properties of isolated single and pairs of substitu- 
tional Mn impurities in GaAs. Apart from helping 
understand the very dilute limit of magnetic semi- 
conductors, these studies are interesting per se^ in 
that individual magnetic impurities in semiconduc- 
tors represent novel nano-magnetic entities with un- 



usual properties and promising applications^ The 
high-resolution STM measurements provide detailed 
information on e.g. the character of the wavefunc- 
tions of single magnetic dopants in semiconductors, 
and the exchange interaction between two nearby 
isolated magnetic impurities interacting via their as- 
sociated itinerant holes. 

Theoretical approaches based on tight-binding 
model a 31 i 32 have so far provided a clear picture of 
many of the STM experimental findings on single 
dopants, including the anisotropic form of the Mn 
acceptor wavefunction and its dependence on the 
Mn spin direction ) 32 ! 33 and the dependence of the ex- 
change coupling between two Mn impurities in GaAs 
on pair orientation and atomic separation! 25 i 31 ' 34 i 35 
Tight-binding calculations are also able to provide a 
description of how the properties of the Mn acceptor 
state depend on the distance of the magnetic impu- 
rities from the surface layer—, which is in qualitative 
agreement with experiment] 29 ! 36 

First-principles calculations based on density 
functional theory (DFT) play an important role in 
the theoretical study of DMSJ^ With the caveat of 
the well-known difficulties of DFT of predicting cor- 
rectly band gaps semiconductor and dealing with 
localized strong correlations at the impurity sites, 
DFT-based first-principle calculations have become 
a very powerful tool to investigate the electronic 
structure and the magnetic properties of DMS. For 
the specific case of the effective exchange coupling 
between two isolated magnetic impurities in bulk 
GaAs, earlier DFT calculation s) 14 ! 15 have provided 
useful information on its microscopic origin and its 
strong anisotropic character with respect to pair ori- 
entation. DFT techniques have also been used to 
simulate STM images of the Mn acceptor states ; 16 ! 37 
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both for interstitial and substitutional impurities, 
yielding results that are qualitatively consistent with 
experiment^ 

A quantity that so far has not been thoroughly 
investigated by first-principle methods is the mag- 
netic anisotropy energy, particularly for the case of 
magnetic impurities located close to the symmetry- 
breaking surface which provides STM access. The 
magnetic anisotropy is a very important quantity, 
particularly when it comes to utilizing these nano- 
magnets in spintronics applications. Indeed, the 
minima in the magnetic anisotropy landscape de- 
termine the direction of the magnetization. In 
DMSs the magnetic anisotropy barriers can be ef- 
ficiently varied with an external electric field, which 
can change the carrier concentration, via the spin- 
orbit interaction (SOI)j£ Therefore it is in principle 
possible to control and manipulate the magnetiza- 
tion direction solely by means of an electric field. 
Tight-binding calculations provide an estimate of 
the magnetic anisotropy for Mn impurities on the 
GaAs surfaced Nevertheless, microscopic calcula- 
tions based on first-principles approach are certainly 
desirable. 

In this paper we present first-principles calcula- 
tions of the magnetic properties for single and pairs 
of substitutional Mn impurities on the (110) surface 
of GaAs, which is the most common cleaved surface 
employed in cross sectional STM studies. We focus, 
in particular on the anisotropic characteristics of im- 
portant magnetic quantities, such as total energy, 
spin and orbital moments and exchange coupling, 
resulting from the interplay of SOI and symmetry- 
breaking surfaces. Since on-site self-interaction cor- 
rections on magnetic impurities are crucial to cor- 
rectly describe the Mn d electronic states in GaAs, 
we carry out our calculations in the framework of the 
Generalized Gradient Approximation + Hubbard U 
(GGA + U) scheme, and we examine the effect of 
on-site correlations on the anisotropic properties. 

We find that the typical magnetic anisotropy en- 
ergy for a single substitutional Mn impurity on the 
(110) plane is of the order one 1 meV, with easy axis 
in the plane. These conclusions are only slightly 
affected by the presence of the Hubbard U term, 
which on the other hand modifies the value of the 
total spin moment. The calculation of the total or- 
bital moment yields very small values, ascribable to 
contributions coming predominately from the accep- 
tor p states and, to a lesser extent, from the Mn d 
states. For two nearby Mn impurities, the magnetic 
anisotropy / atom is of the same order of the single- 
impurity anisotropy. The effective exchange cou- 
pling between the two Mn is typically ferromagnetic 
(i.e., the energetically stable configuration has two 



moments aligned parallel), and strongly anisotropic 
with respect to pair orientation in the surface and 
atom separation. These anisotropics are mainly a 
result of the GaAs crystal structure and symmetry- 
breaking surface. The effect of the SOI on the ex- 
change constant and spin moment is, as expected, 
rather small. 

On the other hand, the exchange coupling is af- 
fected by the Hubbard U term, particularly when 
the Mn atoms are at the shortest possible separa- 
tion. This effect persists at larger separations for a 
Mn pair in the bulk, but it is significantly smaller for 
a Mn pair on the surface, for which the Mn acceptor 
wavefunctions tend to be more localized. 

The paper is organized as follows. In Sec. |TT] wc 
discuss some technical aspects of the DFT calcula- 
tions presented in this paper. Sec. IIII Al contains the 
results of our numerical simulations for individual 
Mn impurities on a (110) GaAs surface. The prop- 
erties on a Mn pair, with a discussion of the effective 
exchange interaction is presented in Sec. IIII Bl Fi- 
nally Sec. IIVI summarizes our work and discusses its 
implications for the physics of magnetic impurities 
in semiconductors and semiconductor spintronics. 



II. COMPUTATIONAL DETAILS 

The density functional theory (DFT) calculations 
in this work are performed using the generalized 
gradient approximation (GGA) with Perdew-Burkc- 
Ernzcrhof (PBE) exchange-correlation functional^ 
For most part of the numerical calculations the 
method of full potential Linearized Augmented 
Plane Wave with local orbitals (LAPW+lo), as im- 
plemented in WIEN2k, is uscd^ Because of the 
time-expensive nature of the plane wave method, we 
have used the SIESTA^ ab-inito package for relax- 
ing the surfaces in our calculations, SIESTA employs 
pseudo potentials and a numerical basis set. The re- 
laxed coordinates are then used as an input for the 
WIEN2k calculations. 

The (110) surface of GaAs for our calculation 
is constructed by cutting the bulk crystal along 

< 110 > direction. The surface supercell consists 
of six layers with 18 atoms at each layer ie total of 
96 atoms. A vacuum of 25 Bohr is added along the 
surface as shown in Fig. [TJ 

To study the effect of an isolated Mn impurity 
on the surface we have replaced one Ga atom from 
the center of the top layer, as shown in Fig. [TJ To 
calculate the anisotropy energy we have calculated 
the effect of SOI for different magnetization direc- 
tions, namely < 100 >,< 001 >, < 011 > and 

< 010 > directions. The anisotropy energy is the 
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FIG. 1: A Relaxed GaAs (110) surface with a Mn impu- 
rity. The arrows indicate different crystallographic axes 
of the surface supercell defined for this work. 



difference between the largest and the smallest en- 
ergies. The spin-orbit coupling in Wien2k is incor- 
porated via a second variational step using scalar 
rclativistic eigcnstates as basis4i 

The strong correlation among the d electrons of 
Mn atoms is accounted for by adding an orbital de- 
pendent potential U for the d shell electrons. The 
value of U for Mn is usually chosen between 3 and 4 
eV to match photoemission spectra and in our cal- 
culation we have used [7 = 4 eV following K. Sato 
et aZ3 

To study the dependence of the exchange energy 
on the pair orientation we replace two Ga atoms by 
two magnetic impurities. One of the impurity atoms 
is kept fixed at the top layer of the supercell (front 
left corner of Fig. [T]) and the position of the second 
one is varied along different crystallographic direc- 
tions, namely [001], [010], [011] and [021]. In this pa- 
per we only consider the case in which the spin mag- 
netic moments of the two impurities are collinear. 
The exchange energy for each pair orientation is cal- 
culated as the difference in energies of the supercell 
when the spins of the impurity atoms are arranged 
parallel and antiparallel. For the anisotropy calcula- 
tion with two impurities we have followed the same 
procedure as we have done for the single impurity 
case. 

Since the size of the surface supercell is large, we 
have used only one k-point for these calculations and 
self-consistency is achieved when total energies are 



converged to within 0.01 mRy or better and charge 
is converged to O.OOle or better. 

III. RESULTS 
A. Single Mn on (110) GaAs surface 

1. Electronic structure 

We start by analyzing the main features of the 
electronic structure when an individual Mn impurity 
replaces a Ga atom on the (110) surface of GaAs. 
A Mn atom introduces a magnetic moment in the 
host material. This can be seen by plotting the par- 
tial Density of States (DOS) for the Mn d-levels as 
shown in Fig. [2] Intra-atomic exchange interactions, 
responsible for Hund's first rule, split the majority 
(up-spin) from the minority (down-spin) d-states, so 
that the former end up - almost entirely - below 
the Fermi level and the latter above, thus giving 
rise to a localized spin moment in the system. In 
Fig. [2] we also show the effect of on-site correlations 
on the d-orbital partial DOS, included via a Hub- 
bard U term in the GGA calculations. Clearly cor- 
relations split majority and minority states further 
apart, pushing the former deeper below the Fermi 
level. As discussed below, correlations also decrease 
the small peak in the partial DOS at the Fermi en- 
ergy as shown in the inset of Fig. [2] 

The second effect caused by a substitutional Mn 
impurity in GaAs is the introduction of a p-like ac- 
ceptor (hole) state in its surroundings. This can be 
seen by plotting the p-orbital partial DOS at the 
Mn nearest-neighbor As sites, as shown in Fig. [3] 
As a comparison, in Fig. G^a) the p-orbital partial 
DOS at an As site is plotted for pure GaAs (i.e., 
when the Mn is replaced by a Ga atom). There we 
can see the expected (full) valence and (empty) con- 
duction band, with a Fermi level in the middle of 
the energy gap that separates them. When the Mn 
is present, the top of the valence band spin polar- 
izes as a result of the hybridization between the As 
p orbitals with the Mn d-orbitals. In Fig. [3f b) we 
can see that the spin- up component (i.e., parallel 
to the Mn spin moment) of the As p-orbital partial 
DOS has a broad peak centered at the Fermi energy. 
The corresponding spin-down component is instead 
pushed below the Fermi energy, causing an effective 
exchange splitting between up and down states of 
the order of ~ 0.6 eV. 

The spin-polarized states just above Ep can be 
identified with the acceptor hole states associated 
with the Mn impurity. These states are clearly the 
result of a moderate hybridization between spin po- 
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larized localized Mn d orbitals with p orbitals pri- 
marily located on the neighbor As atoms. Note that 
the DOS of p character is larger than the d contribu- 
tion. This feature will have an impact on the nature 
of the anisotropy energy discussed below. The inset 
of Fig.f2]shows that on-site correlations decrease the 
majority d partial DOS at the Fermi energy. As a 
result, the hybridization of the hole state with the 
d orbitals decreases, and the acceptor becomes less 
localized on the Mn. 
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FIG. 2: (Color online) The density of Mn d states calcu- 
lated with {7=0 (blue-dark line) and [7=4 eV (red-gray 
line) respectively. The inset shows the DOS near the 
Fermi level. 

It is useful to look more carefully at the effect of 
the Hubbard correlations on the properties of the 
hole state. In Fig. 2] we plot the discrete eigen- 
value energies^ for the states in the surroundings 
the GaAs energy gap, for pure GaAs (red curve) and 
for the case in which a Mn impurity is present (blue 
and green curves). As a result of the Mn presence, 
acceptor states appear in the gap of the host mate- 
rial. Specifically, we find two energy levels located 
a few hundred meV above the GaAs valence-band 
edge. We identify the topmost of these two levels 
as the acceptor state, occupied by a hole. The en- 
ergy immediately below corresponds to the highest 
occupied (by an electron) level. Note that, with the 
exception of one level deep in the valence-band, the 
emergence of these two mid-gap levels is the only im- 
portant qualitative difference that distinguishes the 
spectrum of the system in the presence of the Mn im- 
purity from the one of pure GaAs. As shown in the 
figure, on site U correlations have an effect on these 
two impurity levels, slightly decreasing their energy 
toward the valence band-edge. The sensitivity of 
the acceptor energies to the Hubbard U parameter, 
indicates that these levels are indeed moderately hy- 
bridized with the Mn d orbitals. 
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FIG. 3: (Color online) (Top) The density of As p states 
of pure GaAs surface. (Bottom) The p partial DOS of 
the same As atom when a nearby Ga atom is substituted 
by a Mn atom. The vertical dashed line in both panels 
represents the Fermi energy. The two vertical red lines 
delimit the energy-gap region for pure GaAs, which is 
clearly identifiable in the top panel. To compare the two 
cases, we have aligned the edge of the two conduction 
bands. 
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Eigenstates near the Fermi level 

FIG. 4: (Color online) Eigenvalue energies of eigenstates 
in the surroundings of the GaAs energy gap. The red 
(gray line) is for pure GaAs. The blue (solid black) and 
green (dashed) curves are for the case when a surface Mn 
impurity is present (solid black curve) or not (dashed) 
green. Note the appearance of the hole acceptor state 
in the gap, a few hundred meV above the GaAs valence- 
band edge. The acceptor state energy is moderately af- 
fected by the presence of Hubbard U correlations, a sign 
that this state is in part hybridized with d-states. 
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TABLE I: Square of the wave function, \tpa,tom,ii\ 2 , of dif- 
ferent energy levels near Ef, including the acceptor hole 
state, calculated at the Mn site and its two As nn on 
the (110) surface. The fi — p,d orbital character contri- 
butions are separately specified on the second and third 
column. Here the Hubbard [7 = 0. The first level be- 
longs to the host valence band, while the last level is 
in the conduction band. The bold faced numbers corre- 
spond to the acceptor hole state. The level immediately 
below the hole state corresponds to the highest occupied 
level. Both states are inside the energy gap. See Fig. [4] 



In order to understand further the effect of cor- 
relations on the hole state, in Tables UJ and HI] we 
display the square of the wave function for a few 
levels around the Fermi level at the Mn site and at 
the site of its two As nearest neighbors (nn) on the 
(110) surface. The partial contribution of a given or- 
bital character {p or d) is shown together with their 
sum for a given atom. Table U and [TT] correspond 
to the cases of U = and U = 4 respectively. As 
anticipated above, we can see that the hole state 
comes about from the hybridization of the Mn d or- 
bitals with the p states of its two As nn. We note 
that with U = (see Table H} the hole state is pre- 
dominantly of d character but hybridized with As 
p states. With U turned on (see Table UTJ) , the d 
character of the hole state at the Mn site decreases, 
while the p character weight on the nearby As is es- 
sentially unchanged. This is a clear indication that 
the hole state is pushed out into the interstitial re- 
gion between Mn and its nearby As atoms by on-site 
correlations, namely it acquires a less localized and 
more itinerant character. 

As shown by Mahadevan et al.,— on-site correla- 
tions are known to increase the itinerant character 
of the hole state for Mn impurities in bulk GaAs^. 
To understand the difference in the nature of the 
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TABLE II: Square of the wave function, |V>atom,/i| 
Table HI but with U = 4eV. 



hole state in bulk and surface, we have calculated 
the Mn acceptor hole properties of bulk GaAs. We 
use a 64 atom supcrcell with a single Mn impurity 
in the middle. The occupancies of the hole state at 
the Mn site, its two nn and two next nearest neigh- 
bor (nnn) As sites are tabulated in Table IIIIl It is 
evident that when U is turned on the hole state ex- 
tends over nnn As sites. A similar result is obtained 
by Mahadevan et aZ.— in their bulk GaAs calcula- 
tions. In the surface calculations, although the oc- 
cupancies of the hole state at the nnn As atoms on 
the surface layer slightly increases, the effect is less 
pronounced in the surface compared to the bulk. We 
conclude that for a Mn impurity on the surface, the 
associated acceptor state remains predominantly lo- 
calized in the interstitial region between the Mn and 
the two nn As atoms. Consequently, the hole state 
for a Mn on the (110) is energetically a rather deep 
impurity, compared to the case of Mn in bulk GaAs. 
The localized character of the hole wavefunction for 
a surface Mn and its relatively large binding energy 
are supported by both STM experiment al 29 ' 36 and 
tight-binding calculations.— 

Finally, Tables Q] and [TT] show that on-site corre- 
lations also decrease the Mn e?-orbital contribution 
of the highest occupied level by approximately 50% 
(compare Table HI and . HI)) . As evident from the inset 
of Fig. [2] this effect involves primarily majority spin 
states at the Fermi level. 

The electronic properties the Mn acceptor hole 
state are only weakly affected by the direction of the 
Mn moment direction as a result the SOI. In Fig. [5] 
we plot the difference in the p partial DOS for one of 
the two As nearest-neighbor, for two directions of the 
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FIG. 5: (Color online) Difference in the As p density 
states, for two directions of the magnetization, < 001 > 
and < 010 >, in the presence of SOI and for [7 = 4 eV. 
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TABLE III: Properties of the hole state, as in Tables |H 
Ull but for a Mn in bulk GaAs on the (110) plane. As3 
and As4 are the two next nn As atoms of the Mn on the 
(110) plane. 

magnetization, < 001 > and < 010 > respectively. 
The SOI induced change is only a few percents. We 
will sec below that the magnetic anisotropy energy 
can be explained by this direction dependence. This 
is consistent with tight binding calculations^ for a 
substitutional Mn on the (110) surface, which show 
that the square of acceptor hole wave function is 
anisotropic in space, but does not change signifi- 
cantly with the Mn moment direction. 

2. Spin magnetic moment 

The spin magnetic moment Ms for one Mn on 
the (110) GaAs surface is given in Table ITVl When 
U = the spin moment on the Mn is Ms = 4.05 
in units of the Bohr magneton A close inspec- 
tion shows that there are 4.39 majority-spin and 0.42 
minority-spin d electrons on the Mn atoms. The re- 
maining small contribution 0.08 to the Mn spin mo- 
ment comes from s electrons. The deviation from 
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4.05 


4.07 


u = 


4 


4.46 


4.21 



TABLE IV: Spin magnetic moment on Mn (Ms) and to- 
tal spin moment per unit cell (Mj ot ) in Bohr magnetons. 

the Mn atomic limit of 5 /j,b is due to the accep- 
tor hole being in part localized on the Mn. When 
(7 = 4, the spin on the Mn increases from 4.05 to 4.46 
due to the concomitant increase of the majority-spin 
and decrease of the minority-spin d electrons. The 
increase of the spin moment on the 3d adatom as a 
result of correlations is expected, since the Hubbard 
U term tends to localize the d electrons. Notice how- 
ever that, in this case, the total magnetic moment 
per unit cell Mj ot = 4.21 is smaller that Ms- Indeed 
Mg ot includes also a negative contribution from the 
p orbitals on the nearby As atoms, which are partly 
spin polarized in the opposite direction of the Mn 
d spin moment, as shown in Fig. [3] The results in 
Table IIVI are very weakly affected by SOI. In par- 
ticular, the value of the spin magnetic moment does 
not show any appreciable anisotropy with respect to 
its direction. 



3. Orbital moment 

Next we discuss the orbital magnetic moment. An 
isolated Mn atom has the spin-majority d levels fully 
occupied and the spin-minority states empty. As a 
result the total orbital moment Mr, is zero. Note 
that each of the individual degenerate d orbitals car- 
ries a unit of orbital angular momentum equal to 
mi = 0,±1,±2. When a Mn atom is substituted 
for a Ga atom in GaAs, the crystal field of the host 
material removes the degeneracy of the Mn d lev- 
els and each orbital moment of the perturbed states 
becomes quenched. The effect of SOI is to reduce 
quenching and to restore, at least partially, the indi- 
vidual orbital moments by mixing different d states. 
However, since the spin-majority d manifold is still 
practically fully occupied while the minority-spin is 
empty, we expect that, even in the presence of SOI, 
the total orbital moment Mr, on the Mn atom to be 
very small. On the other hand, an acceptor Mn im- 
purity introduces also an itinerant hole state, which 
is a state with a predominant p character and as such 
should have an orbital magnetic moment. Thus we 
expect that the valence band states, filled up to the 
Fermi level, should generate a net orbital moment 
equal and opposite to the orbital moment of the hole. 
These simple considerations arc fully supported by 
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TABLE V: Orbital d magnetic moment on Mn (Mf) 
and total p orbital moment per unit cell (Mj ot ) in Bohr 
magnetons calculated for two directions of the magneti- 
zation, (001) (easy) and (010) (hard). The minus sign 
indicates that the orbital moment points in the opposite 
direction of the spin orbital moment. 



the DFT calculations. Results are summarized in 
Table El 

We find that the calculated orbital moment of 
the Mn d states is indeed very small and slightly 
anisotropic: for the magnetization directed along 
(001) (which turns out to be the easy direction, see 
next Section) M c [ « 0.015/is; for the magnetization 
directed along (010) (one of the hard axis), is 
about 0.007 /is. The total p orbital moments of the 
unit cell for easy and hard directions are -0.028 and 
-0.016 fiB, respectively These value are larger than 
the d orbital moment but still considerably small. 
The negative moment implies that orbital moment is 
opposite to the spin moment of the system. The or- 
bital moment of the hole state, on the other hand, is 
aligned parallel to the spin moment and, as expected 
form the discussion above, almost exactly cancels 
the orbital moment of the whole cell. The largest 
contribution to the moment comes from the two As 
(14% each) nearest to Mn (5%) and one As in the 
first subsurface layer (12%). For the easy axis, apart 
from these atoms other subsurface layer As atoms 
along [111] direction also have finite contribution to 
the orbital moment, suggesting that the hole state 
is extended along that direction. For the hard axis 
the direction of extension of the hole state remains 
the same but the largest contribution to the p or- 
bital moment comes from the first subsurface layer 
As atom (16.5%), while two surface As atoms con- 
tribute only 7.4% each with negligible contribution 
from Mn. The on-site correlation, U, has very little 
effect on p orbital moment of the whole cell but the 
d orbital moment of Mn becomes vanishingly small 
with U = 4 eV. 

On the basis of these results for the spin and or- 
bital magnetic moment, we can conclude that for a 
Mn impurity on the (110) surface, the total mag- 
netic (spin plus orbital) moment of the system is 
approximately 4 /xb- Assuming that the total mag- 
netic moment is a combined effect of the Mn and the 
acceptor hole, we can imagine to assign an effective 



" spin" J = 2 to this magnetic entity composed of the 
Mn impurity and its hole. The effective spin J can 
be understood as a result of the Mn spin (S = 5/2) 
antifcrromagnctically coupled to the spin of the hole 
state (s = 1/2). For both Mn and acceptor the as- 
sociated orbital moment is small, although not en- 
tirely negligible for the hole. These results can be 
compared with the well-known situation of a substi- 
tutional Mn impurity in bulk GaAs. In that case 
electron-spin resonance^ and infrared spectroscopy 
absorption experiments^ together with theoretical 
considerations, have shown that the the total spin 
of the (Mn + hole) complex is J = 1, and it re- 
sults from the antiferromagnetic coupling of the Mn 
spin moment S = 5/2 and the hole total angular 
momentum j = s + I = 1/2 + 1 = 3/2. Our calcula- 
tions indicate that when a Mn replaces a Ga on the 
Mn surface, the orbital moment of the associated 
acceptor hole is small, most likely due to the fact 
that the high (tetragonal) symmetry experienced by 
the hole in the bulk GaAs is strongly reduced by the 
surface. This reasoning and conclusions are in agree- 
ment with a recent theory^ that identifies the total 
effective " spin" of the (Mn + hole) acceptor magnet 
with a Berry phase Chern number J. Tight-binding 
calculations implementing this theoretical approach 
find that J is 2 for a Mn near the (110) surface, 
due to the quenching of the orbital moment of the 
acceptor hole state. 



4- Magneto-crystalline anisotropy 

We now discuss the magnetic anisotropy energy 
(MAE), defined as the dependence of the ground- 
state energy on the magnetization direction. The 
most important contribution to the MAE is the 
magneto-crystalline anisotropy caused by the SOI. 
For a single substitutional Mn impurity in bulk 
GaAs the tetragonal symmetry of the host lattice 
implies that the MAE is essentially zero. On the 
surface however, the symmetry is broken, and the 
anisotropy is finite. Note that surface relaxation 
and strain (e,g. induced by interfaces with other 
lattices or by external electric fields) can enhance 
the anisotropy considerably. 

Although it is common to calculate the MAE us- 
ing an approximate method known as magnetic force 
theorem^ we have first estimated it by calculating 
the total energy including SOI as a function of dif- 
ferent magnetization directions shown in Fig. [T] The 
anisotropy barrier between two magnetization direc- 
tions, Ms and M' s , is defined as the difference in 
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MAE 


U = 


U = 4 eV 


£((010)) - £((001)) 


1.17 


1.13 


£((100)) - £((001)) 


0.17 


0.16 


£((010)) - £((100)) 


1.00 


0.97 



TABLE VI: Magnetic anisotropy energy (in meV) for 
one substitutional Mn impurity on the (110). 

total energy between these two directions, 

AE A = E(M S ) - E(M' S ) (1) 

In our calculations we find that the total energy 
E(Ms) is minimal when M$ is along the (001) di- 
rection (easy axis) and maximal when Mg is along 
the (010) direction (hard axis), both of which lie in 
the (110) plane. The largest MAE barrier that we 
find from our calculation is 1.17 meV when U = 
and it decreases slightly to 1.13 meV when [7 = 4 
meV. A summary of the magnetic anisotropy energy 
for different magnetization orientations is presented 
in Table ED Notice that the MAE barrier to rotate 
the magnetic moment from the easy axis and align 
it along the < 100 > direction (perpendicular to the 
surface) is only 0.16 meV. 

It is now instructive to look at the MAE form 
the point of view of the magnetic force theorem.— 
This amounts to express the total energy E{Ms) 
as the sum of the magnetization-dependent Kohn- 
Sham band eigenvalues, e„ 

occ 

E(M s ) = Y,£n(Ms) (2) 

where the sum is over the occupied KS eigenval- 
ues. By introducing the (magnetization dependent) 
DOS, we can rewrite E(Ms) as 

E(M S )= [ F (e-e F )p(M s )de, (3) 

where e-Q is the bottom of the valence band and 
ep is the Fermi energy. We can then write the MAE 
barrier as 

MAE = f F (e - e F )[p(M£ ard ) - p{M c s asy )}de . (4) 

J £b 

It turns out that the integrand of Eq. 2] - the 
change in band energy upon varying the direction 
of the magnetization - can be of either sign. When 
these changes are summed up for all occupied states 



we can expect large cancellations, so that what mat- 
ters in the end is the dependence of the band energies 
around the Fermi energy. 

We can now surmise that the magnetic properties 
of the system including the variations of the DOS 
on the magnetization direction, should mainly affect 
the region of the Mn and its immediate surround- 
ings. As we can see in Fig. [5] and Fig. [31 the pro- 
jected DOS at the Fermi energy for the Mn impu- 
rity and its nearest-neighbor As atoms shows that 
the dominant contribution come from the states of 
p character primarily located on the nn As. An in- 
spection of the p DOS for many other atoms in the 
supercell shows that, besides the Mn and its two nn 
As atoms, only the As atoms in the first sub-layers 
closest to the Mn have a finite contribution at Ep, 
which becomes progressively smaller the further the 
As is from the Mn. Thus we can expect that the 
MAE is primarily controlled by the p orbital com- 
ponents of the band energies at Ep located at these 
atoms. 

Indeed, if we calculate the MAE by integrating 
Eq. |4] for all band energies within a few eVs of the 
Fermi level for the Mn and the nearby As atoms, in- 
cluding those As on the sub-surface layers, we obtain 
MAE « 2.8 meV, which differs only by a factor of 
2 from the more precise estimate obtained above. 
The two As atoms nearest to Mn on the surface 
layer contribute about 64% of the MAE, whereas the 
As on the first sub-surface layer and closest to Mn 
contribute about 17%, with Mn d level contributing 
about 18%. 

Furthermore, tight-binding calculations show that 
the MAE obtained by summing up magnetization- 
dependent energy variations for all occupied levels is 
equal and opposite of the anisotropy energy for the 
acceptor hole state only^. If we compute the inte- 
gralMAE holo = (e-e F )[p(M^)-p(M e s asy )]de, 
where ehoic extends up to the limit of the hole-peak 
DOS (see Fig. [3j, we find MAE holc w -2.25 meV, 
which has the expected opposite sign and is very 
close in magnitude to the MAE calculated summing 
all the energy shifts up to Ep. About 68% of the 
hole anisotropy comes from the two As atoms near- 
est to Mn on the surface layer, while the Mn and the 
first sub-surface layer As contributions are 14% and 
16% respectively. 

Note that in the previous section we found that 
a similar property holds for the the orbital moment 
of the system: the total orbital moment, primarily 
located on the Mn, the nn As, and the closest As 
in the immediate subsurface layers, is equal and op- 
posite of the contribution of the hole-state orbital 
moment, located primarily on the same atoms. 

The difference in MAE calculated from the total 
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energy difference (Eq. [T]) and magnetic-force theo- 
rem (Eq. 2]) is expected, since the latter is an ap- 
proximate method, and is due to the fact that DFT 
calculates the groundstate energy of a many-body 
system in a physically meaningful way, whereas in- 
dividual single-particle energies are rather a mathe- 
matical artifact used to solve the Kohn-Sham equa- 
tions. 

In applying the magnetic-force theorem, we have 
disregarded the contribution of all the other atoms in 
the supercell, under the assumption that their con- 
tributions are small around the Fermi level. In fact, 
we have explicitly verified that the p partial DOS 
of these As atoms have a negligible value around 
the Fermi level and the integration of Eq. |4] around 
the Fermi level gives insignificant contribution to the 
MAE. This is, however, not the case if we include 
the band energies deep in the valence band in the 
calculation of the MAE, which appear to give some 
contribution for all atoms. We don't believe that 
this is an indication of a finite contribution to the 
MAE coming from atoms far way from the Mn, but 
rather a sign of the non-sufficient accuracy in the 
computation of the SOI-induced shift in the states 
deep in the valence band. 

At this point it is interesting to investigate more 
closely the connection between the anisotropy prop- 
erties of the total energy and the orbital moment. 
The discussion above already suggests that these two 
quantities should be related. According to the per- 
turbative (in SOI) analysis by P. Bruno^ there is 
an approximate relationship between the MAE and 
the orbital moment anisotropy 

E(M g ) - E(M's) = -| [M L (M S ) - M L (M' S )] (5) 

where £ is the SOI coupling strength. In particular, 
Eq. [5] implies that AEa is proportional to the dif- 
ference in the orbital moment for the spin moment 
in the easy and hard direction. The orbital moment 
calculation clearly shows that the orbital moment 
anisotropy of the hole comes mainly from As p states 
with a small contribution from the Mn d states. 
Using the SOI coupling constant of As, A,4 S =140 
meV^l Bruno formula yields AEf uno =0.35 meV. 
This is approximately a factor of 3 smaller than 
MAE calculated using DFT. The Bruno formula 
has been investigated recently for systems consist- 
ing of magnetic impurities on metal^ and insula- 
tor surfaces.— In some cases the relationship seems 
to work satisfactorily, in other it fails. As far as 
we know a similar analysis for magnetic acceptor 
impurities in semiconductors has not been carried 
out. This situation is different from an ordinary 
magnetic impurity on an insulator for example, in 



that, an important role in the magnetic properties 
of the system is now played by the itinerant accep- 
tor hole. Our analysis shows that orbital moment 
of the hole state is certainly related to the magnetic 
anisotropy energy, as in Eq. [5l albeit with a renor- 
malized SOI coupling strength. This might be an 
indication that our DFT calculations underestimate 
the orbital magnetic moment due to the lack of or- 
bital dependence of the exchange potential^ 

The electronic and magnetic properties of a single 
Mn impurity on the (110) GaAs surface calculated 
here are in general, in good agreement with tight- 
binding based calculations by Strandberg et al^ 
In particular, the value of the magnetic anisotropy 
energy barrier arc pretty close for the two ap- 
proaches « 1 meV. However, details of the magnetic 
anisotropy energy landscape arc different. While in 
our calculation we have both easy and hard axes in 
the (110) plane, Strandberg et al. find an easy axis 
that makes an angle 45° with the surface. The dif- 
ference may stem from the fact that they have used 
much larger clusters and hence their system is much 
more diluted than the system used in this calcula- 
tion. 



B. Mn pair on (110) GaAs surface 

In this section we discuss the properties of a pair 
of substitutional Mn impurities on the (110) surface. 
It turns out that the electronic and magnetic prop- 
erties are quite similar to the ones of individual Mn 
impurities. After a brief summary of these proper- 
ties we will concentrate on discussing the effective 
exchange coupling and the magnetic anisotropy of 
the Mn pair and their dependence on atom separa- 
tion and orientation on the (110) surface. 

1. Electronic structure and spin moment 

The d partial DOS for different pairs of Mn im- 
purities shows a similar pattern as in the case of a 
single impurity, both with U = and U = 4 cV. In 
particular there is rather weak dependence on atom 
separation and pair orientation with respect to the 
surface crystal structure. As for the case of individ- 
ual Mn impurities, the effect of U is to reduce the 
intensity near the Fermi level and to further split the 
main majority and and minority peaks, pushing the 
former below the Fermi level and the latter above. 
Fig.[B]shows the local density of d states of one of the 
two Mn atoms (up spin only) for different pairs on 
the surface, when U =4 cV. While the main majority 
peak is insensitive to the orientation of the Mn pairs, 
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the partial DOS is rather wide at the Fermi level for 
the [010] orientation because of the smaller distance 
between the two Mn atoms for this orientation. 

The properties of the spin magnetic moment (per 
atom) are essentially the same as for individual Mn 
impurities, for all pairs on the surface. For example, 
the Mn spin moment is approximately equal 4 fig 
when U = and increases to 4.5 fig when on-site 
correlations are included by setting U = 4 eV, as in 
the case of a single impurity. 
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FIG. 6: (Color online) Density of d states for different 
orientations of a pair of Mn on (110) surface of GaAs. 
See Fig. [1] for different orientations of the pair relative 
to the crystallographic axes. The inset shows the DOS 
near the Fermi level. 



of the pair and on the pair orientation relative to the 
crystallographic axes of the surface. 

The exchange constant J is largest for the pair 
with the shortest separation along the [010] direc- 
tion, where the two Mn are nearest neighbors, and 
it decays quickly with impurity separation. The 
general trend and order of magnitude of J shown 
in Fig. [7J is in agreement with the results of tight- 
binding calculations, despite these are carried out 
for a much more dilute system, which is supposed 
to model the situation studied in STM experiments. 
Although the smaller size of the surface supercell re- 
stricts our ability to study, in detail, the orientation 
dependence of exchange, it is evident from the fig- 
ure that the exchange constant is larger when the 
pair is oriented along [Oil] direction than when it is 
oriented along [021] direction, even though the dis- 
tance between the Mn atoms is smaller in latter case, 
This is a clear indication of the strongly anisotropic 
nature of the exchange coupling for pairs of substi- 
tutional Mn impurities on the (110) GaAs surface. 
Evidence of this anisotropy was found in recent STM 
experiments^ and is also supported by tight-binding 
calculationsi 25 i 31 ' 34 i 35 We also note that as the dis- 
tance between Mn atoms increases along the [010] 
direction, J becomes negative, i.e. the antiferromag- 
netic coupling of the Mn pair becomes energetically 
favorable. This result implies that a random distri- 
bution of impurities along different directions might 
reduce the overall magnetization of a dilute system. 
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2. Exchange interaction 

A strong exchange interaction between magnetic 
impurities in a semiconductor is essential for its 
successful application as a spintronic device, since 
magnetic properties of semiconductors must be pre- 
served at room temperature and above. In this work 
we have calculated the exchange energy between two 
Mn impurities on (110) surface of GaAs and have 
investigated the effect of correlations and SOI on 
exchange. The exchange energy for each pair orien- 
tation is defined as the difference in energies of the 
supercell when the spins of the two impurity atoms 
are arranged parallel and antiparallel 



J 



(6) 



A positive (negative) J implies that the Mn-Mn 
interactions are ferromagnetic (antiferromagnetic) 
Fig. [Jj shows the dependence of the exchange con- 
stant J on the separation between the two Mn atoms 



-u=o 

-U=4 eV 




1.0 1.2 1.4 1.6 

Impurity separation (lattic constant) 



FIG. 7: (Color online) Dependence of the exchange en- 
ergy constant J on the orientation and Mn-separation for 
a pair of Mn impurities on the (110) surface of GaAs. 

As evident from Fig. the effect of strong cor- 
relations, U, on the exchange anisotropy is rather 
small except for the shortest Mn separation where 
correlations decrease J by approximately 15%. As 
we have seen in the study of one Mn impurity, on- 
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Orientation 


Exchange 


Magnetic 


ol Mn atoms 


energy, J (meV) 


moment 




U=0 


U=4eV 


U=0 


U=4eV 


[010] i 


267.9 


223.7 






[001] 


33.3 


35.5 






[021] 


3.2 


8.9 


4.0 


4.4 


[010] 2 


-12.4 


-7.6 






[Oil] 


19.4 


19.8 







TABLE VII: Exchange energy and magnetic moment for 
a Mn pair on the (110) GaAs surface. The subscripts " 1" 
and "2" in the [010] direction indicate respectively the 
nearest and 2nd nearest Mn atoms along that direction. 



site correlations have a small tendency in dclocaliz- 
ing the acceptor hole around the Mn. According to 
theory, the acceptor hole is expected to mediate the 
exchange interaction among the magnetic moments 
of Mn pairs. Except for the shortest impurity sepa- 
ration, the enhanced itinerant character induced by 
correlations is still of limited range in space to have 
a noticeable effect on the exchange constant. Never- 
theless, a close inspection of the numerical values of 
J (see Table lLUB2[) shows that, with the exception 
of the pair with the shortest separation, the effect of 
U is always to increase J This is consistent with 
the fact that exchange mediating-holc becomes more 
dclocalized as a result of on-site correlations. 

Although the results shown in Fig. [7] are obtained 
in the absence of SOI, we have checked that SOI 
has little influence on J. For instance for the [010] 2 
orientation J increases by 1-2 meV depending on 
the direction of magnetization, for both U = and 
U = 4 eV. 



largest J, is also the one with the smallest MAE. For 
all other Mn pairs for which J fluctuates around a 
considerably smaller value, the corresponding MAE 
fluctuates around 2 meV. It is remarkable that ex- 
actly the same behavior is found in tight-binding 
calculations^ where the value of the anisotropy for 
a Mn pair is also in quantitative agreement with the 
results presented here. 

In the presence of correlations U = 4 eV, the 
MAE, in general, decreases except for the case when 
for the pair with the shortest Mn separation. This 
behavior is again consistent with the fact that, with 
the exception of the shortest pair separation, on-site 
correlations slightly increase the exchange constant 
J. 
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FIG. 8: (Color online) Surface anisotropy for a pair of 
Mn impurity on (110) surface of GaAs. 



3. Magnetic anisotropy 

The MAE barrier for a pair of Mn atom on GaAs 
(110) surface varies from 0.44 meV to 2.67 meV de- 
pending on the orientation of the pair relative to 
the crystallographic axes, as shown in Fig.[5J For all 
pairs except for the pairs along the [010] orientation, 
the easy axis is found to be along < 001 > and the 
hard axis is along < 010 > direction. For the pair 
oriented along [010], while for the nn the situation 
is inverted (i.e., the easy axis is along < 010 > and 
the hard direction is along < 001 >), for the next 
neighbor pair the easy axis is still along < 001 > but 
the hard axis changes to the < 100 > direction. 

Our calculations also show an interesting correla- 
tion between the MAE and the exchange constant 
J. Comparing Figs. [7] and [8] we note that the Mn 
pair with the shortest Mn separation, which has the 



IV. CONCLUSIONS 

In this work we have investigated the electronic 
and magnetic properties of substitutional Mn im- 
purities on the (110) surface of GaAs, using den- 
sity functional theory implementing a full potential 
LAPW+lo method. 

A key point of our analysis has been the investiga- 
tion of of the Mn acceptor hole state, which appears 
in the GaAs energy gap. The properties of the hole 
state allow us to understand the salient features of 
a single Mn impurity in GaAs. As for a bulk Mn 
impurity, when the Mn is on the surface, the associ- 
ated hole state is strongly anisotropic in space, but 
much more localized around the Mn and its two As 
nn than in bulk. 

On site electron correlations on the Mn, included 
within a GGA + U framework, reduce the hybridiza- 
tion of the hole state with the Mn d orbitals, and 
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render the character of the wave function relatively 
more p-like and therefore slightly more delocalized. 
The spin magnetic moment of the Mn atom is ap- 
proximately 4/iR, when U = 0, consistent with a hole 
partly localized on the impurity. The Mn d moment 
increases to 4.5/xb when U = 4 eV, as a result of 
the more delocalized nature of the hole state. How- 
ever the spin moment of the whole unit cell remains 
closer to 4/ib- 

SOI induces a small dependence of the p par- 
tial DOS at Ep on the direction of the magneti- 
zation, primarily for the As atoms closest to the 
Mn. With some caveats (see discussion at the end 
of Sec. IIII A 4[) . this dependence accounts for a cal- 
culated magnetic anisotropy energy (MAE) of ap- 
proximately 1 meV and can be closely related to the 
anisotropy of the hole state. 

The presence of SOI yields a small orbital mag- 
netic moment in the opposite direction of the spin 
moment, which can again be associated to the or- 
bital moment of the hole state. The dependence of 
the As p partial DOS on magnetization directions in- 
duces a very weak anisotropy in the orbital moment. 
We have seen that the Bruno formula relating the 
MAE to the orbital anisotropy is qualitatively satis- 
fied, but underestimates the MAE. 

When two Mn impurities are substituted for two 
Ga atoms on the (110) surface we find that the most 
energetically stable configuration is typically the one 
where the two Mn magnetic moments are oriented 
parallel to each other. The exchange energy J be- 
tween two Mn impurities on the surface strongly 
depends on the orientation of the pair relative to 
crystallographic axes. On the other hand, the SOI- 
induccd dependence J on the direction of the mag- 
netization is very small. The effect of on-site cor- 
relations on J is largest for the pair along the [010] 
direction with nn Mn, where the exchange is reduced 
by 15%. For all the other pairs, correlations always 
enhance exchange, although the effect is typically a 
few percents. Finally, the MAE of a ferromagnetic 
Mn pair on the GaAs surface is slightly larger than 
that of a single impurity, except for the shortest pair 
for which the MAE is 0.2 meV/Mn. With U added, 
the MAE decreases in general. 

The results presented here are in good agreement 
with recent tight-binding calculations carried out for 
individual and pairs of Mn substitutional impurities 



on the (110) surface of GaAsj 32 ' 35 ' 44 In particular, 
the MAE of a single impurity is in both cases on 
the order of 1 meV. Similarly the exchange energy 
for a Mn pair displays the same rapid decay with 
Mn separation and anisotropic behavior with pair 
orientation. 

The orbital magnetic moment of the acceptor hole 
is very small, as a result of the surface that breaks 
the tetragonal symmetry of bulk GaAs, again con- 
sistent with tight-binding calculations 4i This sug- 
gests that for a Mn impurity on the surface the "to- 
tal spin" of the magnetic center composed of the 
Mn core and the associated spin-polarized accep- 
tor hole should be S'xot = 5/2 — 1/2 = 2, result- 
ing from antiferromagnetic coupling between the Mn 
spin (Smii = 5/2) and the hole spin (sh)- If this 
"spin" is subject to an anisotropy landscape with 
energy barriers of 1 meV, it would take a magnetic 
field on the order of 10 T to revert its direction. Al- 
though for a Mn on surface the acceptor wavefunc- 
tion is rather insensitive to the direction of the Mn 
magnetic moment, this does not seem to be the case 
for impurities located in the first 10 sub-layers below 
the surface. In particular, the acceptor partial DOS 
on (110) surface accessible from STM experiments is 
predicted to depend strongly on the direction of the 
magnetic moment. This would imply the possibility 
of controlling the STM tunneling current by rotating 
the Mn acceptor with a magnetic field. So far exper- 
iments carried out with magnetic fields up to 7 T do 
not find any dependence of the conductance for dif- 
ferent directions of the magnetic fields If the MAE 
barrier of the impurity is indeed 1 meV or larger, 
such fields might simply be not strong enough to 
rotate the direction of the spin. 
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